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The spin-half XXZ model on the linear chain and the square lattice are examined with the ex- 
tended coupled cluster method (ECCM) of quantum many-body theory. We are able to describe 
both the Ising-Heisenberg phase and the XF-Heisenberg phase, starting from known wave functions 
in the Ising limit and at the phase transition point between the XF-Heisenberg and ferromagnetic 
phases, respectively, and by systematically incorporating correlations on top of them. The ECCM 
OS \ yields good numerical results via a diagrammatic approach, which makes the numerical implemen- 

OS ■ tation of higher-order truncation schemes feasible. In particular, the best non-extrapolated coupled 

OS ■ cluster result for the sublattice magnetization is obtained, which indicates the employment of an im- 

proved wave function. Furthermore, the ECCM finds the expected qualitatively different behaviours 
£^ ■ of the linear chain and the square lattice cases. 
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I. INTRODUCTION 



The extended coupled cluster method (ECCM) has not previously been applied to lattice spin systems, unlike 
the normal coupled cluster method (NCCM) [0-13 1, which is a restricted version of the ECCM at a given level of 
approximation, and which has been widely implemented for these systems. The primary aim of this paper is to 
apply the ECCM to the spin-half anisotropic Heisenberg (or XXZ) model, in order to obtain numerical results for 

os ■ 

the ground-state energy and the sublattice magnetization, and thereby to investigate the usefulness of the method 
in the study of quantum phase transitions. The ECCM and NCCM are ah initio techniques of microscopic quantum 
many-body theory, generically known as the coupled cluster method (CCM). The ECCM, by contrast to the NCCM, 
completely characterizes a system in terms of a set of basic amplitudes, all of which are linked-cluster quantities. 

The two-dimensional spin-half XXZ model is expected to have a second-order phase transition at the isotropic 
point, A = 1, where the system is in a unique critical phase. The transition is expected to be accompanied by some 
change of symmetry of the ground-state wave function. In particular, Laughlin jl4j speculates that the physics of 
the isotropic point can be understood as a gauge theory with massless excitations. In general, theoretical study of 
' the isotropic point requires some prior knowledge of the phases on either side of the isotropic point. In practice, the 
ordering and symmetry of the known wave functions in these phases influence the predicted state at the isotropic 
point. Therefore, the second aim of this paper is to examine how the CCM, in particular, is affected by this universal 
problem in the study of phase transitions. 

The general XXZ model Hamiltonian has Z(2) ® U(l) xy symmetry, except at the isotropic point where it has SU(2) 
symmetry. For the square-lattice case the ground-state wave function is expected to yield broken U(l) xy symmetry 
in the XF-like region, and there is broken Z(2) symmetry in the Ising limit. Approaching the isotropic point from 
either side with the CCM yields different ground-state wave functions at the isotropic point itself, which is due to the 
special nature of the isotropic point. 

Any CCM calculation on spin systems involves adding correlations between spins, on top of those already contained 
in a separately chosen model or reference state, in order to produce the true quantum-mechanical ground-state wave 
function. Therefore, the important physical characteristics of a system are incorporated into the CCM by the choice 
of the model state and by the inclusion of particular configurations for the correlations. Hence, we can now refine and 
reformulate our second aim to be an examination of the effect of the choice of model state. In particular, we shall be 
interested in the interplay between the symmetry of the model state and the symmetry of the Hamiltonian influencing 
the ground-state wave function. This is exemplified by the CCM treatment of the linear-chain case, which yields 
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artificially broken U(l) xy and SU(2) symmetry. This can only be due to the choice of model state, since exact results 
show that the symmetry of the XXZ model Hamiltonian is not broken in the XF-like region and at the isotropic point. 

The primary aim of this paper is motivated by numerical evidence, from a number of techniques, of long-range 
order (LRO) at the isotropic point for the square-lattice case. Numerical results from techniques such as spin-wave 
theory [p~5|— p~T|| , high-accuracy quantum Monte Carlo (QMC) simulations M], and series expansions yield a 

sublattice magnetization at the isotropic point of approximately 61-62% of the classical value arising from perfect 
Neel ordering. However, it is still possible that each of these techinques may underestimate the quantum fluctuations 
in a similar fashion, since they all start from an ordered state 

Rigorous results from the Bethe ansatz |2^-Q for the spin-half anisotropic Heisenberg antiferromagnet on the 
linear chain provide a measure of the effectiveness of the ECCM. However, we know in advance that the nature of the 
transition in the linear chain is very subtle, and quantum fluctuations present in this case are known to destroy Neel 
LRO completely at the isotropic point p5| . Therefore, a priori, one would expect the ECCM to be more effective for 
the square-lattice case, where the ground state will be closer to the classical uncorrelated state. 

Numerical results for the ground states of models in one and two dimensions are expected to be qualitatively 
differe nt jP6j . By contrast with the NCCM, which has previously been applied, for example, to the spin- half XXZ 
model |27|] and the non-linear sigma model [E8[, the ECCM is expected to yield numerical evidence of this difference. 
As similar calculations are performed for any lattice dimensionality in the CCM, the qualitatively different behaviour 
of the solutions that we report in this paper for the two cases using the ECCM is not simply an artefact of the 
technique. 

In a diagrammatic implementation of both CCM techniques, the ECCM, at any level of truncation, produces 



diagrams in greater abundance and of greater complexity than the NCCM, as shown in Sec. IV, Unlike in its 
NCCM counterpart, it is possible within the ECCM formalism to define a spin-spin correlation function which is 
fully consistent with the corresponding definition of the magnetic order parameter, at all orders of truncation. Such 
properties of the ECCM make it attractive for describing correlation effects in physical systems, and particularly to 
study their quantum phase transitions and quantum order. 

The CCM has been used very successfully |j|-[l3) to calculate the zero-temperature properties of a wide variety 
of extended many-body systems, including, for example, atoms and molecules, nuclear matter and finite nuclei, the 
electron gas, lattice gauge field theory, as well as spin-lattice systems [^7],|9]-33 . In this paper, we present evidence 
from numerical results that demonstrates the particular superiority of the ECCM over the NCCM, in practice, to 
study such global properties of spin systems as their quantum (zero-temperature) phase transitions. 



II. EXTENDED COUPLED CLUSTER METHOD FORMALISM 



Since detailed descriptions of the CCM formulation have been given elsewhere [|l|-|13| , only the essential components 
will be given here. Hubbard |3^ | was one of the first to emphasize the importance of an exponential parametrization 
of the exact ground-state ket wave function \ip) of an interacting many-body system, 

|y,)=e s |0) ; S = J2^iC\ , (2.1) 
in terms of a model stat e \4>) which is not orthogonal to the exact ground-state ket wave function The correlation 



operator S in Eq. (2.1) is decomposed solely in terms of a complete set of mutually commuting, linked, multicon- 
figurational creation operators C\ defined in terms of a complete set of many-body configurations {/}. These are, 
in turn, defined by a set-index /, which, for the spin-lattice systems under consideration, describes the set of spins 
which are flipped with respect to those contained in the suitably chosen, normalised, model state \<p). The prime on 
the sum in equation ( |2.l| ) excludes the null set, 1^0, corresponding to the identity operator, Cq = 1. We note that 
Ci\(f>) = 0, MI ^ 0, and hence the state obeys the intermediate normalisation condition, (0 1 -0) = 1. 

Arponen 0] has shown how two distinct CCM parametrizations of the exact ground bra state {ip\ can be given, 
which yield the NCCM and the ECCM respectively. The NCCM parametrization of the ground bra state is given as: 

(TP\ = (<f>\Se- s ; S = l+Y,hC! , (2.2) 

whereas the corresponding ECCM parametrization is given as: 

= (<f>\e s " e~ s ; S" = J2siCj . (2.3) 
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Both S and S are constructed wholly in terms of multiconfigurational destruction operators, defined with respect to 
the model state \<j>). These are simply the hermitian adjoints, Cj, of the corresponding creation operators, Cj, in Eq. 



(2.1). Furthermore, both parametrizations satisfy the normalization condition, (ip\ip) = 1- Although (ip\ = / 
formally, this relation may not be preserved when truncations are made, as explained below, in cither parametrization. 
The energy functional He in the ECCM formulation is thus given by: 

H E = {^\e s " e- s He s \(f>) . (2.4) 
The similarity transformed Hamiltonian e~ s He s may be expressed as the usual nested commutator expansion, 

e~ s He s = H + [H, S] + 1 {[H, S] , S] + ■ ■ ■ . (2.5) 



Since S is composed wholly of mutually commuting creation operators, the similarity transform in Eq. (2.5) only 
retains terms in which all amplitudes si are linked to the Hamiltonian. Furthermore, provided that the Hamiltonian 
is finite order in the single-body operators, the expansion always terminates at finite order. Therefore, once the 
correlation operator S is approximated, no further truncations are necessary. 

Formally, the bra-state parametrizations in both the ECCM and NCCM formulations still preserve hermiticity. 
Explicitly, we have the relations, 



® S = M* =W^W) • (2 - 6) 

which are a consequence of satisfying the normalisation condition (ip\ip) — 1. Although the manifest hermiticity is 
usually sacrifi ced at any level of truncation, the distinct bra-state parametrizations produce fully linked expectation 



values in Eq. (2.4), as they incorporate the similarity transform, and both S and S are composed wholly of destruction 
operators. 



The double exponential structure of the ECCM formalism implies that Eq. (2.4) can be expressed in terms of 
the amplitudes {s/,Sj} via a double similarity transform, and that both of these sets of ECCM amplitudes are 
linked-cluster quantities. However, the NCCM formalism only allows a single exponential structure in terms of the 
amplitudes {sj}. Although the amplitudes {s/} are linked-cluster quantities, the amplitudes {s/} are not. 

An arbitrary expectation value A can thus be expressed in the ECCM and the NCCM respectively as: 



A E ee (<f,\e ff 'e- 8 Ae 8 \4>)=A B (i ai }, {sty , (2.7) 
A N ee (<j>\Se- s Ae s \<j>) = A N ({ Sl },{Sj}) . (2.8) 
In practice, either version of the CCM can only be implemented by truncating the expansions in Eqs. (2.1)-(pT3j) 



by retaining only a finite or infinite subset of the complete set of configurations indices {/}. In this paper the SUB2-n 
approximation scheme is employed, which retains all configurations with up to two-body correlations between spins 
which are no more than a certain distance apart, specified by the index n. The details of the approximation scheme 



are given explicity in Sec. IV where the diagrammatic approach is described. 



All ground-state properties can be determined in the ECCM and the NCCM respectively, by the stationary principle: 

dH E dH E dH N _ _ 3H N 

-r. — — u — „ , —z — — u — , (z.y) 

osj dsj osi osj 

in the case where |^) and {tp\ correspond to the ground state. We remark, however, that due to the lack of manifest 
hermiticity between {ip\ and at a given level of truncation, the resulting stationary values of He and Hn are not 
necessarily upper bounds to the ground-state energy. 

We note that the NCCM equations at a given order can formally be extracted from the corresponding ECCM 

equations at the same order by truncating the expansion of the exponentiated correlation operator e s at first order 
and performing the substitution Sj — > Sj. In practice, the resulting ECCM equations from a particular truncation 
scheme are highly non-linear and, therefore, of greater complexity than their NCCM counterparts. However, the 
diagrammatic representation of the formalism which we present here makes tractable the numerical implementation 
of high-order approximation schemes. 

Very importantly, observables in the ECCM, which quantify the global behaviour of a system in terms of its 
long-range order, obey the cluster property. Thus, for example, we have the very general relation, 
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lim A r B T , = A r B r , , (2.10) 

|r— r | — >oo 

where A r and B T i are single-body operators acting at denned sites r and r , respectively. This condition is preserved 
by the ECCM due to the exact multiplicative separability of both the bra-state and ket-state parametrizations in 
the corresponding large-distance limit, however the index set {/} is truncated. Consequently the long-range order of 
physical systems can be examined unambiguously via the ECCM parametrization of their correlation functions. As 



the NCCM bra-state parametrization is not multiplicatively separable, however, the physical condition in Eq. (2.1C) 
will not in general hold at arbitrary levels of truncation within this scheme. 



III. TRANSFORMATION OF THE SPIN-HALF XXZ MODEL HAMILTONIAN FOR THE LINEAR 

CHAIN AND SQUARE LATTICE 

A. The CCM approach to the spin-half XXZ model Hamiltonian 

The spin-half XXZ model Hamiltonian is given by: 

ff = |^[A^ + ^ + ^J , (3.1) 

where af,a = x,y, z, are the Pauli spin matrices, A is the anisotropy parameter, and the summation is over all the N 
lattice sites denoted by i and over each of the nearest-neighbour vectors denoted by p. We note that the linear chain 
and the square lattice are both bipartite lattices which can be split into two identical sublattices, which we denote as 
the A- and -B-sublattices. Thus, each nearest-neighbour site to a site on the A-sublattice is on the B-sublattice, and 
vice versa. 



In the Ising limit, A — > oo, the A-dependent term in the Hamiltonian of Eq. (3.1) becomes dominant and the 
classical z-aligned Neel state is the eigenstate that yields the lowest energy. The z-aligned Neel state has nearest- 
neighbour spins ordered antiparallel to one another in the z-direction, as illustrated below: 

\<t>) = ®|T>* 011)/ , (3.2) 

keA leB 

in a notation in which the z-axis points vertica lly upwards. However, for all finite values of A the terms in the 



x- and y-directions in the Hamiltonian in Eq. (3.1) come into play, and the z-aligned Neel state is no longer an 
eigenstate of the Hamiltonian. The ground state for even N now consists of a particular linear combination of all 
possible configurations with N/2 up-pointing spins and N/2 down-pointing spins. All the configurations apart from 
the classical state which are present in the exact state are considered to be quantum fluctuations upon that state. At 



the isotropic point the Hamiltonian in Eq. (3.1) becomes rotationally invariant, such that the expectation value of 
an arbitrary spin is the same for any direction. 

For — 1 < A < 1 the true classical ground state can be any one of an infinite number of degenerate Neel states with 
nearest- neighbour spins restricted to align antiparallel to one another in any direction in the XY plane. Schematically, 
we write: 

10} = ®l<->* (g)H>i , (3-3) 

keA leB 

choosing, say, the x-direction to be the alignment axis, and in a notation in which the positive x-axis points horizontally 



to the right. We note that there exists a trivial transformation of the Hamiltonian in Eq. (3.1) at the phase transition 
point, A = —1, between the XF-Heisenberg and ferromagnetic phases (henceforth referred to as the ferromagnetic 
point), with the cc-aligned Neel state as eigenstate, to the ferromagnetic Heisenberg Hamiltonian, with the x-aligned 
ferromagnetic state as eigenstate. Therefore, at A = — 1 the x-aligned Neel state is the true ground state and this 
point is chosen to be the initial point of the CCM calculation in the entire regime — 1 < A < 1. One imagines that 
the x- aligned Neel model state is close to the true ground state in the neighbourhood of this point. 
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B. Transformation of the Hamiltonian in terms of a canted model state 



In order to produce a ground-state solution that can describe both the Xy-likc (—1 < A < 1) and Ising-like (A > 1) 
regions of the XXZ model, a canted model state is introduced. The canted model state, \<p c ) = |0(£a,£b)}, consists of 
spins on the ^4-sublattice pointing in one particular direction and spins on the _B-sublattice also pointing in another 
particular direction, defined by the spin- half spin vectors £a and £b, respectively. 

In order to calculate the expectation values of arbitrary observables in the ECCM, an arbitrary rotation of the local 
spin axes is performed about the y-axis on each sublattice resulting in a notional rotation of the spins in the canted 
model state to the down position in the direction of the negative z-axis, 

N 

UmU,^B)) = (g)\i)i = \F) , (3.4) 

i=i 

where U is a product of unitary matrices, which causes the spins in the canted model \4> c ) to undergo a passive 
rotation such that th ey a ll point downwards in the rotated local frames. Hence, the ECCM expectation value for an 



observable A in Eq. (2/7) can be expressed with respect to the unrotated canted model state in the form, 

A = (4)^11^ e s " e~ s (UAU^) e s U\4> c ) . (3.5) 
The rotation matrices for the canted model state on the A-sublattice and B-sublattice are given by: 



Uj = exp \~i6j-j = cos 11 - * sin y-f j ° v ; J = A B , (3.6) 

where each rotation is chosen to be performed about the y-axis in the XZ-planc. This involves no loss of g ener ality, 
since all directions in the XF-plane are equivalent, due to the U{l) xy symmetry of the Hamiltonian in Eq. (3.1). 



Now, the transformed Hamiltonian arising from Eq. (3.1) can be written as 



H T ee UHU^ = +\ ]T [A (u A alu\) (u B afU B ) + (lJ A a%U^) (u B afU B 

(k,i) 

+ (u A a y k U A ^j (u B afU^] , k e A, I e B , (3.7) 

where k and I are nearest neighbour vectors. 

There are two degrees of freedom, 8 A and 9b, present in the transformed Hamiltonian H T , which are conveniently 
expressed via an equivalent set of relative and total orientation parameters, 

a={9 B -9 A ) ■ f3 = -{9 A + 8 B ) . (3.8) 

Particular values of a and (3 lead to specific canted model states, as shown in Fig. 1. The single-spin creation and 
destruction operators are defined as cr ± = \ (o x ± icr v ) for all sites on the lattice, once the spins have all been rotated 
into the down position, as given in Eq. (3.4). Thus, H T in Eq. (FO) can be expressed in the form: 



where p, q £ {z, +, — } and T vq are functions of a, (3 and A. Due to the hermiticity of the Hamiltonian the factors T pq 
in equation ( |3.9| ) satisfy the relation, T +z = T* z . Moreover, the spatial symmetry with respect to which lattice sites 
the spin operators act upon, yields the relations: 

T z+ = T +z ; T,_ = T_ 2 ; T+_ = T_+ . (3.10) 

Finally, the transformed Hamiltonian with a canted model state can be expressed as: 

rT = +1 E { \ K 1 + A ) cos a + ( J - A ) cos P - 2 1 ( a trt +P + °I°i +P ) 



i,p 



~ [(1 + A) cos a + (A - 1) cos/3] o\a z i+p 
- [(1 + A) cos a + (1 - A) cos (3 + 2] cr+ar +p 

-(l-A)sin/3(a.f CT + p + ^- p )} . (3.11) 
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C. A mean-field calculation 



A mean-field calculation provides a guideline for the model-state analysis in the correlated CCM calculation, which 



will be performed in Sec. |VI 
of the canted model state 



We restrict our mean-field calculation here to finding an optimal state of the form 
(£a>£b))j m which no other correlations are present apart from those between the two 



sublattices implied by the angles a and j3. Consequently, the mean-field energy E, which is given by: 



E 



\H\<j> c ) = (F\H T \F) 



yields the following expression: 



— = ^[(A + l)cosa + (A-l)cos/3] 



(3.12) 



(3.13) 



where N is the number 
Hamiltonian in Eq 
T 

x zz 



of lattice sites and z is the co-ordination 



(3.9) yields the only non-zero, diagonal contribution from Eq. (3.12) 
with respect to the angles a and (3 yield the solutions for the energy E, which are shown in Fig. 
shows that the model states which make up the ground state in the mean-field case are the z-aligned Neel state for 
A > 1, the x-aligned (or y-aligned) Neel state for — 1 < A < 1, and the z-aligned ferromagnetic state for A < — 1. 



number of the lattice. The factor T zz in the 
The stationary values of 
Figure □ also 



IV. THE DIAGRAMMATIC APPROACH 



A. The Diagrammatic Technique 

The terms arising from a CCM calculation have been represented diagrammatically in previous work J2^j35|. For 



example, Roger and Hetherington |29| represent the terms arising from the NCCM energy functional diagrammatically 
simply for convenience. Unlike our representation, shown in Fig. ||, their diagrams do not include the destruction 
operators from the bra-state parametrization. Hence, diagrams are retained in Ref. pS[ ], which do not contribute to the 
energy functional, as shown by the NCCM diagrams in Fig. || On the other hand, Harris uses a representation 
where the algebra for performing the CCM with respect to spin states is reformulated in terms of particle-hole 
states. This particular representation was employed for the purpose of identifying cancelling terms and to aid the 
simplification and systematization of the algebra. 

By contrast with these earlier diagrammatic approaches, the originality of the diagrammatic approach p8| used 
in the present work appears at the level of the numerical calculations. Since, the diagrammatic representation 
requires less computer memory and CPU time than the algebraic representation, higher-order truncation schemes 
can be implemented. There are several other advantages of using this particular approach: applying the formalism is 
straightforward; the terms which contribute the most in any calculation can readily be pinpointed; and the correlation 
function and sublattice magnetization can be expressed as sets of the diagrams which appear in the CCM energy 
functional arising from H T . 

Our diagrammatic technique involves casting the combinations of creation and destruction operators which appear in 
the CCM in a diagrammatic form, as shown in Fig. ^. The CCM SUBn approximation scheme retains all configurations 
up to and including n-body correlations. Therefore, the SUB2 scheme retains all the one-bod y cor relation terms, 



denoted by Si in Fig. 3, that are allowed to arise from the terms in the Hamiltonian H T in Eq. ( 3.11 ) which contain 
a single-body creation or destruction operator, as well as all the two-body correlations, denoted by 52 in Fig^. [| The 
SUB1 scheme retains only the one-body correlations. In the present work, the correlation operators S and S defined 



in Eqs. ( [2.l| ) and (2.3) are truncated via the SUB2 scheme as shown in Fig. |3J. 

As the equations that arise from the SUB2 scheme are not analytically soluble, a further restriction is placed on 
the maximum range of the two-body correlations from the SUB2 scheme. Therefore, the partial SUB2 approximation 
known as the SUB2-n scheme is implemented. In the current diagrammatic formulation the SUB2-n scheme retains 
only those two-body correlation coefficients b r , where r runs over the n distinct vectors within a reference box of given 
size. The vectors arc defined to be distinct after allowing for all lattice symmetries. For example, the linear-chain 
case has the property that all correlation coefficients arising from vectors of the same length are identical. However, 
for the square-lattice case, vectors of the same length which are not otherwise identical under the lattice symmetries 
yield distinct correlation coefficients. For example the vectors (5,0) and (3,4) have the same length but are distinct. 

The one-body terms in Si from the SUB1 scheme simply transform a particular model state <fi to an improved 
model state \(f> ), where \<jj) is in the neighbourhood of \cf> ), via the transformation \(f> ) = e Sl \<fi). The one-body terms 
from the SUB1 scheme are rigorously zero when the chosen model state is the same as the ground-state mean-field 



G 



states, as the SUB1 scheme corresponds to a mean-field calculation. Thus, the particular cases of the z-aligned and 
x-aligned Neel model states for A > 1 an d — 1 < A < 1 respectively, as shown in Fig. 0, from the Hamiltonian 
H T in the local rotated axes in Eq. (3.11), yield rigorously zero one-body coefficients. More general canted model 
states yield very small one-body contributions. Therefore, in order to keep the number of diagrams under control, we 
restrict the one-body terms to give, at most, a linear contribution, which corresponds to the NCCM SUB1 scheme. 
Hence, the restricted SUB2 bra-state parametrization can now be expressed as: 

(4.1) 




exp 



'i+r 



(4.2) 



Once the energy functional (2.4) arisi ng from H T is calculated, the ground-state properties can be derived by 
applying the stationary principle of Eq. (2J5). An example of how the diagrammatic calculation is performed can be 
show n by taking the expectation value of the term of crf +p , which is a constituent of H T . The nested commutator in 
Eq. (2J3) yields the following expression: 



1 



0"iCT* +p + [ofof +p , S]+- [[ofof+p, S] , S] 



(4.3) 



The correlation operator, 5*2 in Fig. ||, which is described as a ket line, flips two spins in the model state \F) to 
the up position. The nested commutator in Eq. (4.3) only allows non-zero contributions, if the ket line is connected 
to the term of of +p , in the Hamiltonian. Therefore, at least one of the spins is flipped on either of the lattice sites i 
and i + p. Consequently, a spin is flipped on one other site of the lattice. 

The bra lines arising from e only yield non-zero contributions with the state described above if the final state is 
one of all down-pointing spins. As the bra lines flip up-spins to down-pointing spins, the ends of the bra lines must 
be connected to ket lines or to a term in H T that contains creation operators. Otherwise, the destruction operators 
at either end of the bra line annihilate the state. The combinatorics of both the bra lines and the ket lines in devising 
the diagrams is taken care of by the CCM formalism. The ket lines arising from the nested commutator and the bra 

lines arising from e s have counting factors that remove the possibility of obtaining diagrams which are equivalent 
under lattice symmetries. 



B. The diagrams arising from the ECCM formulation for macroscopic quantities 

1. The diagrams arising from the energy functional 

The diagrammatic technique applied to the ECCM energy functional arising from the transformed Hamiltonian 
H T yields the diagrams shown in Fig. |]. 

2. The sublattice magnetization 



The sublattice magnetization, or the order parameter in the Ising-Heisenberg and XY"-Heisenberg phases, respec- 
tively, is defined by, 



keA 



keA 



keA 



N 



A > 1 



-1 < A < 1 



(4.4) 



keA 



The inclusion of the minus sign ensures that M z and M x are positive in our rotated Neel basis, where U z and U x are 
the unitary operators, which rotate the local spin axes of the z-aligned and cc-aligned Neel model states respectively 
to give the state \F). The sublattice magnetization for A > 1 is solely that of the z-aligned spins M z , and for 
— 1 < A < 1 the sublattice magnetization is that of the x-aligned spins M x in the present model state formulation. 
The formal reason for this behaviour of the sublattice magnetization is given in Sec. [v|. 
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As the forms of the sublattice magnetizations M z and M x for the respective regions shown in Eq (4.4) are equivalent, 
they arise from the same diagrams, which are shown in Fig. |^. Local quantities in the CCM formulation yield size- 
extensive results. The diagrams for the ECCM sublattice magnetization are the same as those for the NCCM. Even 
though they are local quantities, the order parameters arise from summing correlations over the whole lattice, and in 
this way they are able to describe the ordering of individual spins throughout the system. 



3. The correlation function 

In order to study critical phenomena it is important to study a particular correlation function which is consistent 
with the definition of the order parameter. The spin-spin correlation function g^ is defined by: 

So'V) = Wk a k+r) > a,P = x,y,z , (4.5) 
in terms of which we define the correlation function g a '^(r) as follows, 

rt^rt-wxO • ( 4 - 6 ) 



where r is restricted such that k and k + r are both on the A-sublattice. The translational invariance of the lattice 

J3 

i-\-r 



implies, if a = f3 on a particular sublattice, (erf) = {&i +r )- Hence, the correlation function satisfies the fundamental 



property: 

iflS"» = flE<^>! = (M a ) 2 , (4.7) 



r — >oo 



where M a is the order parameter and i is summed over the sublattice of N/2 lattice sites. The expression for go(r) 
in terms of the rotated canted model state is given by: 

go(r) = (UaM+rtf) + {U^^rf) + {Ua z k ot +r U^) , (4.8) 



where U is the unitary operator in Eq. ( p.q ). Clearly the form of the expression in Eq. Q4.8| ) is similar to the form of 
the corresponding energy functional. The same diagrams arise for go(r) as for the energy functional, which is shown 
in Fig. ^. However, the links can now be of arbitrary length, r. The property in Eq. (4.7), implies that in the 



l imit r — > oo, surviving diagrams from go(r) can be interpreted as being topologically disconnected. Conversely, Eq. 



(4.6) implies that the correlation function g(r) consists of the topologically connected diagrams arising from go(r). 
Examining the topology of diagrams arising from the functional provides a straightforward way to deduce the terms 
that yield the sublattice magnetization and the correlation function in an ECCM calculation. 

The relevant correlation functions for the regions which encompass the Ising-Heisenberg and XY"-Heisenbcrg phases 
respectively are given by: 

g>*(r) = (U'vlai^U*)-(U'o%U*)(U'ai+ r U't) ; A > 1 , 

9 ™(r) = (U x aM +r U^)-{U x a%U^)(U x a* k+r U^) ; -1 < A < 1 . (4.9) 



In the limit r — > oo, the correlations between spins die off as shown in Eq. (4.7). The spin-spin correlation function for 
the region -1 < A < oo has the form ±(of of +r ) as (U z afa z +r U z ^) — (U x (j x (j x +r U x ^) = ±(of of +r ), where the sign 



is positive if both lattice sites are on the same sublattice and negative otherwise. Therefore, the correlation function 

\ +P ) 



consists of diagrams similar to the connected diagrams aris ing from the term (crfah_ p ) in the energy functional H T 



The diagrams arising from the correlation function in Eq. (4^) are shown in Fig. |6| 

The connected diagrams which appear in the NCCM correlation function are equivalent to those in Fig. || with the 
exception of those indicated in the box. However, there are also disconnected diagrams in the NCCM formulation 
of the correlation function. This occurs because one of the diagrams that arises from the square of the sublattice 
magnetization in Fig. ^| is disconnected and corresponds to a higher-order term that cannot appear in the spin-spin 
correlation function, due to the NCCM bra-state parametrization. Consequent ly, t he cancellation of the disconnected 



diagram cannot occur when the correlation function is calculated in Eq. (4.9). However, the ECCM bra-state 
parametrization yields higher-order terms than the NCCM, which allows this cancellation to take place. Hence, the 
correlation function is solely composed of connected diagrams, which possesses the fundamental property of Eq. (|4.7| ). 
In this way the importance of the cluster property, which is intrinsic in the ECCM formalism in Eq. ( 2.10| ), can be 
seen for the correlation function. 
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V. THE SYMMETRIES OF THE MEAN-FIELD MODEL STATE 



Criteria for a good choice of model state in terms of the maximum overlap between the model state and the ket 
state have been derived by Kummel |36| . One of these criteria specifies that there should be no contribution from the 
one-body terms (or the SUB1 scheme) in the CCM calculation. The model states arising from the ground-state mean- 
field calculation satisfy this criteria. From the canted model state, the mean-field calculation yields the ground state 
for the ^-aligned Neel model state in the region A > 1 and the ^-aligned Neel model state in the region — 1 < A < 1, 
as shown in Fig. ^. However, it can be shown, by applying a theorem first enunciated by Xian p7[ to these mean-field 
model states, that the ket state \ip) cannot exhibit symmetry breaking at A = 1. 

Xian's theorem [^| states that the CCM equations provide at least one solution which guarantees the symmetry of 
the model state, if this symmetry is one of those belonging to the Hamiltonian. The theorem is proved by employing 



a symmetry operator A, which has the properties that the Hamiltonian in Eq. commutes with A and that the 
model state is an eigenstate of A. 

In the case of the z-aligned Neel model state, the symmetry operator A z is given by the rotation operator, 

JV 



A* 



mn(f)<7?] , (5.1) 



where i runs over all lattice sites. The order parameter for an arbitrary direction in the XY-plane is defined by 
M xy = M x cos a + M v sin a. The action of the rotation operator A z on the operators whose expectation values on 
the sublattice yield M xy is given by: 

A z \a x cos a + o v sma)A z = a x cos (a + x) + & y sin (a + x) ■ (5.2) 

The only solution that guarantees the required invariance of M xy under rotation is M x = and M y = 0, such that 
the U(l) xy symmetry of the z-aligned Neel model state is preserved. Consequently, the theorem |37j shows that \ip) 
cannot exhibit broken U(l) xy symmetry. 

Similarly the discrete symmetry operator for the particular case of the x-aligned Neel model state is given by: 

N 

V (g) of • (5.3) 

i=l 

The transformation of the order parameter M z to — M z with A x implies that M z = 0. Thus, the Z(2) symmetry of 
the x-aligned Neel model state is preserved and consequently \ip) cannot exhibit broken Z{2) symmetry. Hence, it 
follows that the z-aligned Neel model state only describes the region A > 1, where the ground state does not have 
broken U(l) xy symmetry. Similarly, the x-aligned Neel model state only describes the region — 1 < A < 1, where the 
ground state does not have broken Z(2) symmetry. 



VI. NUMERICAL RESULTS 

The diagrams in the ECCM energy functional corresponding to the Hamiltonian H T , which arise from the SUB2-n 
approximation scheme in Fig. |^, are used to obtain the numerical results for both the linear-chain and square-lattice 
cases. The qualitative behaviour of the numerical results differs for the two cases. This difference can only arise from 
the different multiplicities of the diagrams, which themselves are due to the symmetries of the lines in the diagrams 
differing in these cases. 

From the implementation of the NCCM SUB1 scheme it can be seen that the coefficients k and k provide a 
negligible contribution to the overall calculation, where they are not rigorously zero for the x-aligned and z-aligned 
Neel model states. Therefore, with the benefit of hindsight, implementing the ECCM SUB1, which only introduces 
additional higher-order terms, would not be much different to the NCCM SUB1. 



A. The spin-half XXZ model on the linear chain 



There are both exact results for the case of the linear chain and more appropriate numerical techniques than 
the CCM, such as the density matrix renormalization group and spin-wave theory, for examining more general one- 
dimensional spin chains. However, the presence of exact results makes the linear chain an important case to gain 
insight into the behaviour of the numerical results from the approximated ECCM formalism. 
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The mean-field calculation shows that both the z-aligned and x-aligned Neel model states introduce complete 
z-aligned and x-aligned Neel ordering respectively, at the isotropic point. Figures and || show that the ECCM 
SUB2-n scheme fails to pick up all the quantum fluctuations in the physical system that result in the LRO being 
completely destroyed at the isotropic point. However, the exact order parameter from the linear-chain case exhibits 
subtle behaviour, which, in principle, would require the implementation of extremely high-order ECCM approximation 
schemes. 

The ECCM numerical solutions, for the quantities M x and E g , from truncation schemes of orders higher than 
SUB2-2 with the x-aligned Neel model state, at the initial point A = —1, do not exist for A > —1, as shown in Figs. 
| and g. This occurs for any scheme higher than SUB2-2 because there is a degeneracy in the numerical solutions 
at A = —1, which makes it a poor choice of starting point. An alternative choice of initial point for the CCM 
calculation is the isotropic point, since the rotated Hamiltonians obtained from both the x-aligned and z-aligned Neel 
model states are equivalent at this point. These solutions turn back and become unphysical. The points at which the 
solutions become unphysical are known as the terminating points, and are shown in Figs. |^ and ^. 

The terminating points from the SUB2-n scheme based on the s-aligned Neel model state, at the orders that 
have so far been obtained, converge weakly to 1.025 ± 0.005 as n — > oo, which can be interpreted as the SUB2 
terminating point. This result was obtained from the best possible mean square fit of 1/n ' 563 . It indicates that the 
x-aligned Neel model state from the SUB2 approximation scheme would not yield a solution in the XF-like region, 
where it is considered to be a better model state than the z-aligned Neel model-state. This is surprising since higher 
orders include more of the long-range correlations in the system, and consequently should reveal more of the physical 
behaviour of the system. The exact ground state of the linear chain in the XF-like region possesses U(l) xy symmetry, 
which indicates that the x-aligned Neel model state, as the sole input in the technique, breaks this symmetry causing 
solutions to terminate in the region A < 1. However, the numerical results from the ECCM for the ground-state 
energy and sublattice magnetization M x are good at low orders. 

Figure [l(] shows that the low-order ECCM SUB2-n solutions have the same qualitative behaviour as the exact 
ground-state energy solution, where a maximum deviation in the energy occurs at the isotropic point. Numerical 
results from the NCCM LSUBn scheme [j32|j yields a numerically accurate ground-state energy at the isotropic point 
from an extrapolation, as shown in Table. |l The LSUBn scheme retains only those configurations in the correlation 
operator S which contain any number of spin flips with respect to the model state over a localized region of n 
contiguous sites, and which are compatible with the restriction Sj~ = | 2i=i °f = 0- We note that this restriction 
follows from the facts that Sf. commutes with H and that the ground state is expected to be in the Sj, = sector. 

B. The spin- half XXZ model on the square lattice 

There are no exact results for the square-lattice case, which therefore presents a significant challenge for numerical 
techniques, particularly for the study of the nature of the expected transition at the isotropic point. Unlike the 
linear-chain case, the square-lattice case is believed to possess LRO at the isotropic point. Evidence for this comes 
from a number of techniques as shown in Table |J. The quantum fluctuations in the square-lattice case should be 
smaller than those in the linear-chain case, due to the higher co-ordination number of the lattice. For this reason one 
expects that the CCM should be better at deducing the physical behaviour of the square-lattice case. 

The numerical solutions for M z and E g , from truncation schemes at orders higher than SUB2-12 with the z- 
aligned Neel model state, terminate before the isotropic point, as shown in Figs. |TI| and [l2| respectively. The SUB2 
terminating point of the CCM calculations based on the z-aligned Neel model state can be accurately determined by an 
extrapolation, where the best mean square fit of 1/n 2 gives a SUB2 terminating point at Ak 1.03903±0. 00077. Since 
the z-aligned Neel model state only breaks the discrete Z(2) symmetry of the ground state, the SUB2-n terminating 
points exhibit this good convergence. 

Similarly, solutions starting at A = — 1 from the SUB2-20 scheme with the x-aligned Neel model state for E g and 
M x terminate in the XY"-likc region, as shown in Figs, [l^ and |TJ respectively. The SUB2-20 approximation scheme is 
the highest we have so far investigated, as well as being the only scheme to yield a terminating solution in the XY"-like 
region. By contrast, at lower orders terminating points occur for A > 1. Therefore, although the SUB2 terminating 
point cannot be accurately determined, it is expected to be closer to A = — 1 than the SUB2-20 terminating point 
at 0.78. Consequently, it is not possible to obtain a solution from high-order SUB2-n schemes that can describe the 
neighbourhood of the isotropic point on the A < 1 side. 

Unlike in the linear-chain case the existence of a solution from high-order schemes with the ^-aligned Neel model 
state in the square-lattice case indicates that the continuous U (l) xy symmetry of the ground state is broken. However, 
the bad convergence of the terminating points from the x-aligned Neel model state in the square-lattice case may 
indicate that the ground state is almost symmetric and, hence, a model state with U(l) xy symmetry would be a more 
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appropriate starting point for a CCM calculation. 

The SUB2-12 scheme is the highest order SUB2-n scheme that yields the sublattice magnetization M z at the 
isotropic point, because at higher orders the solutions terminate before the isotropic point. The SUB2-12 scheme yields 
a value of M z k, 0.689 of the saturation value. This is an improvement upon the NCCM SUB2 value of M z w 0.81, 
and the best NCCM result from the LSUB8 scheme J32j, which yields a value of M z « 0.705, when compared to other 
numerical techniques, as shown in Table. ||. The extrapolated numerical value of the magnetization from LSUBoo, 
M z « 0.622, is the result of a quadratic fit in 1/n with a restricted set of three data points. The best linear fit in 1/n 
with the same data points yields M z w 0.646 ± 0.002 

Quantitatively, the CCM yields accurate values for the ground-state energy E g compared to other numerical tech- 



niques as shown in Table. The NCCM SUB2 scheme |27j at the isotropic point yields a ground-state energy of 
Eg « —0.651, which the ECCM SUB2 scheme improves upon, via an extrapolation with values that converge well, 
with a numerical result of E g « —0.667. The calculation of the ground-state energy from the NCCM LSUBrt schemes 
p2| are extremely accurate when compared to the best Monte Carlo result |l8| where the limit n — > oo is equivalent 
to including the complete set of many-body configurations. 



C. The canted model state in the correlated ECCM SUB2-n calculation 



The effect of the canted model state, which has two degrees of freedom, in a CCM calculation is examined here. 
In practice, the spins in the canted model state are restricted to point in one direction on a particular sublattice. 
Nevertheless, for the spin-half XX Z model, the CCM with the canted model state provides a good analysis of the 
effect of the model state in these calculations. In terms of the numerical calculations, performing the CCM with the 
canted model state is a new approach, since the earlier CCM calculations for spin systems have chosen pre-determined 
model states with no free parameters, and with the choice based on either prior knowledge or classical behaviour. 
The use of the canted model state allows the extra freedom to partially tailor the model stat e to the specific value 
of the anisotropy parameter. This is achieved by the application of the stationary equation (2J5) and the variation 
with respect to the angles a and (3, which define the model state, where the minimum energy is automatically chosen 
for a and (3 at each value of the anisotropy A. Therefore, we employ this particular form for the numerical CCM 
calculation, in the hope that it will chose model states which lie closest to the true ground state. 



1. The behaviour of the canted model state with varying angles 



For the calculations with varying anisotropy, the angles a and (3 are free to vary. The model states that arise from 
the ECCM SUB2-n calculation in this case are the same as those from the mean-field calculation, as shown in Fig. |^, 
for — 1 < A < oo. This is also illustrated by Fig. |lj in the neighbourhood of the isotropic point. As discussed in Sec. 
0, these mean-field model states yield the CCM ground state that preserves either the U(l) xy or Z(2) symmetry of 
the Hamiltonian. 

Kiimmel's work |3(J provides a guideline to examine the effect of the mean-field model states in the CCM calculation, 
if translational invariance of the rotated model state \F) is assumed. It can be shown that the overlap condition for 
the mean-field model states is maximized if the sum, over the whole lattice, of the two-body correlations, which is 
defined in Fig. || as S2, obey ^2 r b r < 1/2. This result is extended to show that J2 r b r — * 1/2 indicates a phase 
transition point. As A — » 1 this condition holds for the z-aligned Neel model state in the linear-chain case and only 
approximately holds for the square-lattice case. For the square-lattice case, this may indicate the numerical instability 
of the solution as A — > 1 and consequently the need for a model state with greater overlap in this region. 

The solution for the energy at the isotropic point has a continuous degeneracy in terms of the angle /?, for a fixed 
value of a = it which yields the set of all model states with Neel symmetry. However, away from the isotropic point 
there is a unique minimum in the solution. Therefore, in practice, by choosing a fixed angle /?, we break the SU (2) 
symmetry at the isotropic point. In doing so, LRO is introduced in the direction specified by the angle f3, which is 
diminished by the correlations from a CCM calculation. 



2. The behaviour of the canted model state with fixed angles 

The effect of particular model states that require the inclusion of one-body terms in the CCM calculation are 
examined here. In the Ising-like region, for the angle (3 set to 7r, the angle a is fixed to various values which define 
model states in the neighbourhood of the z-aligned Neel state, i.e., a ~ 7r. As a is decreased, each numerical solution 
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for the energy in the square-lattice case, from the highest-order ECCM SUB2-n scheme, is successively higher and 
terminates earlier than the previous solution. Furthermore, for model states defined by a < 7r/2, a solution cannot 
be recovered. 

These non-mean-field model states, in the neighbourhood of the z-aligned Neel model state, require one-body 
terms in the CCM calculation. Hence, the ground state can exhibit broken U{l) xy symmetry |37|. However, these 
model states are bad choices in comparison to the symmetry-preserving mean-field model states, due to the poor 
qualitative and quantative behaviour of the numerical solutions for the energy. Moreover, due to the termination of 
these numerical solutions well before A = 1, broken U(l) xy symmetry does not physically occur in this region. Thus, 
the symmetry-breaking nature of these poor choices of model state cannot be utilised. It seems that non-mean-field 
model states with one-body correlations try to approach the mean-field model states and, in doing so, do a worse job 
at describing the system. 

3. The non-trivial behaviour of the canted model state in the Ising limit 

We now examine the effect of the canted model states in the Ising limit. As shown in Fig. |l4|, for the Ising-like 
region of the linear-chain case the lowest energy, and therefore the optimum model state to yield the ground-state 
energy, is the z-aligned Neel state. Initially, in the Ising-like region the periodicity of the solution for the ground-state 
energy in terms of the angle (3 is 2ir, when a is fixed to 7r, as shown in Fig. [lq for the linear-chain case, and with 
minima at ir mod(27r). In the Ising limit the solution asymptotically develops a periodicity of it in the angle (3. The 
states at (3 = mod(27r), which define the x-aligned Neel state, become degenerate with those at (3 — ir mod(27r), 
which define the z-aligned Neel state. Interestingly, this analysis reveals that in the Ising limit, the z-aligned Neel 
ground state is obtained from the the x-aligned Neel model state with large one-body coefficients k. 

VII. DISCUSSION AND CONCLUSIONS 

The diagrammatic technique has made possible the implementation of high-order SUB2-n schemes in the ECCM 
formulation, since it requires less computer memory and CPU time than performing the CCM with algebraic methods. 
This is because the multiplicties of the diagrams are only determined once and can thereafter be used throughout the 
calculation. For the case of the square lattice, accurate numerical results in comparison to other numerical techniques 
have been obtained for the region — 1 < A < oo. In particular, at the isotropic point the most accurate value for the 
order parameter of M z » 0.689 for any non-extrapolated CCM calculation has been obtained. This suggests that 
the ECCM SUB2-n scheme describes the wave function better than previous NCCM formulations at similar levels of 
truncation. 

The qualitative behaviour of the numerical results indicate that the present CCM treatment requires further im- 
provement for the study of phase transitions. In the case of the linear chain the ECCM fails to describe the system for 
|A| < 1. For the square-lattice case there is a failure to describe the region near to the isotropic point. Alternatively 
the manner in which the numerical solutions become unphysical provides further knowledge of the system and the 
ECCM. 

It is intuitively clear that a good choice of model state in a CCM calculation should be as close as possible to to 
the exact ground state or, alternatively, share the same underlying symmetries as the exact ground state, depending 
on the system under consideration. The high-order SUB2-n ECCM numerical results show that the important facet 
of a good choice of model state is dependent on the extent to which the system is ordered at the isotropic point and 
in the XY"-like region. The linear-chain case with both the z- and ^-aligned Neel model states, and the square-lattice 
case with the re-aligned Neel model state, indicate that when either LRO does not exist or the value of the order 
parameter may be quite small, it is important to chose a model state which does not artificially break the symmetry 
of the Hamiltonian. 

Alternatively, the square-lattice case with the z-aligned Neel model state suggests that if the system has broken 
symmetry, the large quantum fluctuations that result from the massless excitations could also cause solutions to 
terminate. In this case, it is more important to chose a model state which is as close as possible to the exact ground 
state. This problem of correlations becoming too large seems to occur generally in quantum many-body techniques. 
Advantageously, introducing broken symmetry into such a system yields the benefit of not having to perform high- 
order calculations to approximate the true ground state well. 

As shown above, our analysis, as it stands, requires prior knowledge of the system. Although other numerical 
techniques |l^-|20[| suggest that the square-lattice Heisenberg model possesses LRO, the value of the order param- 
eter might be small. Therefore, the best CCM approach for this system is to define a model state with a flexible 
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parametrization, so that it does not artificially break the symmetry of the Hamiltonian, and in so doing can give 
indications of possible symmetry breaking by examining the effect of correlations on the order parameter. In order 
to achieve this end the model-state parametrization should yield a mean-field solution in which the order parameter 
smoothly and continuously tends to zero at the isotropic point. In this way the ECCM formulation can be extended 
over the whole regime to include the long-range correlations at the critical point. 
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TABLE I. The ground-state energy per spin for the one-dimensional chain at the isotropic point, A = 1, under various CCM 
schemes, compared with the exact results. 



NCCM SUB2 


ECCM SUB2 


NCCM LSUBoo 


Exact 


-0.419 


-0.433 


-0.443149 


-0.443147 



TABLE II. The sublattice magnetization M z for the two-dimensional square lattice at the isotropic point, A = 1, as a 
fraction of the saturation value, from various CCM approximation schemes, a series expansion, a Monte Carlo simulation, and 
spin- wave theory (SWT). 



ECCM SUB2-12 


NCCM LSUB8 a 


NCCM LSUBoo". 


series expansion 13 


Monte Carlo c 


SWT d 


0.689 


0.705 


0.622 


0.614 ±0.002 


0.615 ±0.005 


0.606 



a Reference 
b Reference 
c Reference 
d Reference 



TABLE III. The ground-state energy per spin for the two-dimensional square lattice at the isotropic point, A = 1, from 
various CCM schemes, compared with the results of a Monte Carlo simulation. 



ECCM SUB2 


NCCM LSUBoo a 


Monte Carlo 10 


-0.667 


-0.66968 ± 0.00004 


-0.66934 ± 0.00003 



a Reference 
b Reference 



2 

18] 



FIG. 1. Various model states canted in the XZ-plane, defined by the angles a and j3, are shown in the original global 
co-ordinate frame, in a notation in which the positive a;-axis points to the right and the positive z-axis points upwards. Note 
that the states with a — ir are all Neel states and those with a = are all ferromagnetic states. 

FIG. 2. The stationary energy eigenvalues, E/N, as a function of the anisotropy, A, for the mean-field case. The ground-state 
energy is denoted by a thick line. 

FIG. 3. The diagrammatic representation of the ECCM formalism. The dashed lines denote the interaction terms from the 
Hamiltonian. The circle and the the cross denote respectively a creation operator and a destruction operator. The straight 
solid lines are known as the ket lines and the wiggly solid lines are known as the bra lines. 

FIG. 4. The diagrams arising from the ECCM energy functional (2.4) using the Hamiltonian H T , which yields the 
ground-state energy. The factors denoted by the letters are composed of functions of the angles a and /3 and the anisotropy 
A. The factor z denotes the co-ordination number of the lattice. The diagrams that appear up to the solid line are those from 
the NCCM formulation. 



FIG. 5. The diagrams constituting the sublattice magnetizations M z and M x 



FIG. 6. The diagrams for the correlation functions, g zz and g* 
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FIG. 7. The sublattice magnetization in the z direction, M z , as a function of the anisotropy A for the linear-chain case, 
from the ECCM SUB2-n schemes using the 2-aligned Neel model state, a mean-field calculation, and the exact results. 

FIG. 8. The sublattice magnetization in the a;-direction, M x , as a function of the anisotropy A for the linear-chain case, 
from the ECCM SUB2-n schemes and a mean-field calculation with the a;-aligned Neel model state. For truncation schemes 
higher than the SUB2-2, solutions turn back, so only the physical solution is shown with a solid circle indicating the turning 
back point, also known as the terminating point. It is known that there is no LRO over the whole region, — 1 < A < 1. 

FIG. 9. The ground-state energy per spin, E g /N, as a function of the anisotropy A for the linear-chain case from high-order 
ECCM SUB2-n schemes using the x-aligned Neel model state, compared with the exact results. 

FIG. 10. The ground-state energy per spin, E g /N, as a function of the anisotropy A for the linear-chain case, from the 
low-order ECCM SUB2-2 scheme using the z-aligned and rr-aligned Neel model states in the regions A > 1 and — 1 < A < 1 
respectively, and the exact results 

FIG. 11. The sublattice magnetization in the 2-direction, M z , as a function of the anisotropy A for the square-lattice case, 
from the high-order ECCM SUB2-n schemes using the z-aligned Neel model state. The results converge at high orders of the 
SUB2-n truncation scheme, thus giving the full SUB2 result, which corresponds to n — > oo. 



FIG. 12. The ground-state energy per spin, E g /N, as a function of the anisotropy A, from the 2-aligned and x-aligned Neel 
model states in the regions A > 1 and — 1 < A < 1 respectively, for the square- lattice case. The terminating points of the 
solutions from the ECCM full SUB2 and SUB2-20 approximation schemes, in the regions A > 1 and — 1 < A < 1 respectively, 
are indicated by solid circles. 

FIG. 13. The sublattice magnetization in the x-direction, M x , as a function of the anisotropy A for the square-lattice case, 
from a mean-field calculation and high-order ECCM SUB2-n schemes using the a;-aligned Neel model state. 

FIG. 14. The ground-state energy per spin, E g /N , as a function of the orientation of nearest-neighbour spins, (3, in a model 
state, which is restricted to the canted model state space, for the linear-chain case, when the relative angle a is it (Neel state). 
The examination of the model-state space is performed in both the J^F-like and Ising-like regions in the neighbourhood of the 
isotropic point, using the converged ECCM full SUB2 values. 

FIG. 15. The ground-state energy per spin, E g /N, as a function of the overall orientation of nearest-neighbour spins /3, when 
a is fixed at it, for the linear-chain case. The examination of the model-state space is performed in the Ising-like region for 
increasing values of the anisotropy, using the converged ECCM full SUB2 values. 
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